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Abstract 

We consider the dynamics of a harmonic crystal in d dimensions with n compo- 
nents, d, n arbitrary, d,n > 1 , and study the distribution fit of the solution at time 
t £ R. The initial measure fio has a translation-invariant correlation matrix, zero 
mean, and finite mean energy density It also satisfies a Rosenblatt- resp. Ibragimov- 
Linnik type mixing condition. The main result is the convergence of fit to a Gaussian 
measure as t — > 00 . The proof is based on the long time asymptotics of the Green's 
function and on Bernstein's "room-corridors" method. 
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1 Introduction 



Despite considerable efforts, the convergence to equilibrium for a mechanical system has re- 
mained as an extremely difficult problem. It has been recognized early on that for an infinitely 
extended system, possibly on top of local hyperbolicity, the flow of statistical information to 
infinity serves as a mechanism for relaxation. The two prime examples are the ideal gas and 
the harmonic crystal. We consider here the latter case. In the harmonic approximation the 
crystal is characterized by the displacement field u(x) , where x G V , V is a regular lattice 
in H d , and u(x) G R n with n depending on the number of atoms in the unit cell. The 
field u(x) is governed by a discrete wave equation. We will consider arbitrary d,n and for 
notational simplicity set T = Z d . 

Our motivation to return to a well studied model is to a much wider class of initial 
measures than before. This project requires novel mathematical techniques. They have 
been developed for the wave and Klein-Gordon equation on H d in [6] - [8], but the discrete 
structure poses extra difficulties. 

Let us briefly comment on previous work. In [14] a general criterion is given which ensures 
mixing and Bernoulliness of the corresponding mechanical flow. Thereby the convergence to 
equilibrium is established for initial measures which are absolutely continuous with respect 
to the canonical Gaussian measure. In [14] moments of the displacement field are studied. 
This allows to reduce the spectral analysis of the Liouvillean flow to the spectral properties 
of the dynamical group defined on solutions of finite energy. Since the crystal is assumed to 
be homogeneous, these spectral properties are determined by the dispersion relations u>k(9) , 
k = 1, n . The Liouvillean flow is mixing and even Bernoulli, if, except for crossing points, 
each u>k{0) is a real- analytic function which is not identically constant. In particular, the 
Lebesgue measure of the set {9 G T d : \/uk(9) = 0} is equal to zero. In [20], for the case 
d — n — 1 , initial measures are considered which have distinct temperatures to the left and 
to the right. In [2], again d — n — 1 , the convergence to equilibrium is proved for a more 
general class of initial measures characterized by a mixing condition of Rosenblatt- resp. 
Ibragimov-Linnik type and which are asymptotically translation-invariant to the left and to 
the right. 

The detailed stationary phase analysis of [2] does not directly generalize to d > 2 . Rather 
we have to develop a novel 'cutoff strategy' which more carefully exploits the mixing condition 
in Fourier space. This approach allows us to all d with in essence the same conditions for the 
dispersion relations as in [14]. Our extension requires the technique of holomorphic functions 
of several complex variables. 

In parenthesis we remark that, for the ideal gas, R.Dobrushin and Yu.Suhov [3] first 
realized the importance of a mixing condition on the initial measure. In [9] it is replaced by 
the condition of finite entropy per unit volume thus establishing convergence whenever the 
specific particle number, energy, and entropy are finite. No such general result seems to be 
available for the harmonic crystal. 

We outline our main result and strategy of proof. The displacement field u(x) is the 
deviation of the configuration of crystal atoms from their equilibrium positions. Assuming 
them to be small and expanding the forces to linear order yields the discrete linear wave 
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equation, 

u(x,t) = -^2 zd V(x - y)u(y,t); u\ t =o = u (x), u\ t =o = v (x), xeZ d . (1.1) 

Here u(x, t) = (ui(x, t), . . . , u n (x, t)),u = (u i, . . . , u 0n ) e R n and correspondingly for v . 
V(x) is the interaction (or force) matrix, (Vki(x)^j, k,l = l,...,n. The dynamics (1.1) is 
invariant under lattice translations. 

Let us denote by Y(t) = (Y°(t),Y\t)) = (u(-,t),u{;t)) , Y = (YgX) = MO, *><>(•)) • 
Then (1.1) takes the form of an evolution equation, 

Y(t) = AY (t), t e R; Y{0) = Y . (1.2) 

Formally, this is the Hamiltonian system since 

AY = j(l ° i y = JVH(Y), ■/=(_° 1 J)- (1-3) 

Here V is a convolution operator with the matrix kernel V and if is the Hamiltonian 
functional 

H(Y) = ±(v,v) + ^(Vu,u), Y = (u,v), (1.4) 
1 1 

where the kinetic energy is given by -(v,v) = -V] vd \v(x)\ 2 and the potential energy 

by -(Vu,u) = ~^2 x y( - zd (V(x — y) u {y)i u { x )) j (' ) ') being the real scalar product in the 
Euclidean space R n . 

We assume that the initial datum Y is a random element of the Hilbert space H a of real 
sequences, see Definition 2.1 below. Y is distributed according to the probability measure 
/i of mean zero and satisfying the conditions S1-S3 below. Given t G R, denote by /i t the 
probability measure for Y(t) , the solution to (1.2) with random initial data Y . We study 
the asymptotics of fit as t — > ±oo . 

The correlation matrices of the initial data are supposed to be translation- invariant, i.e. 
for i,j = 0, 1 , 

Q%(x,y):=E(Y*(x)®Y ( i(y))=qi ) \x-y), x,yeZ d , (1.5) 

though our methods require in fact much weaker conditions. We also assume that the initial 
mean "energy" density is finite, 

e := E[\u (x)\ 2 + \v (x)\ 2 } = tr g 00 (0) + trg n (0) < oo, xeZ d . (1.6) 

Finally, it is assumed that the measure fio satisfies a mixing condition of a Rosenblatt- resp. 
Ibragimov-Linnik type, which means that 

Yq(x) and Y Q (y) are asymptotically independent as \x — y\ — > oo. (1.7) 

Our main result is the (weak) convergence of the measures /i t on the Hilbert space Ti a with 
a < —d/2 , 

—r Hoc as t — > oo. (1.8) 
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//oo is a Gaussian measure on H a . A similar convergence result holds for t — > — oo . Explicit 
formulas for the correlation functions of the limit measure //oo are given in (2.18) - (2.22). As 
an application of the results, we show that the initial "white noise"- correlations provide the 
limit measure fi^ which coincides with the Gibbs canonical measure with the temperature 
~ e • Respectively, /i^ is close to the canonical measure if the initial correlations are close 
to the white noise. 

To prove the convergence (1.8) we follow general strategy [2, 4, 6, 7]. There are three 
steps: 

I. The family of measures // t , t > , is weakly compact in H a with a < —d/2 . 

II. The correlation functions converge to a limit, for i, j = 0, 1 

Q?(x,y) = jY i (x)^Y^y) f , t (dY)^QZ(x,y) as t - oo. (1.9) 

III. The characteristic functionals converge to a Gaussian one, 

:= J exp(i(Y, *)MdY) - exp^Q^, *)} as t - oo. (1.10) 

Here * = $ 1 )gD = D®D, D = C (Z d ) <g> R™ , where C (Z d ) denotes the space of 
the real sequences with finite support, (Y,^) = J2i=o,iJ2xez d (Y l ( x ) > ^(xfj and Qoo is the 
quadratic form with the matrix kernel (Q^x, y))ij=o,i , 

eoo(*,*)= E E (QUx,y),nx)®* j (y))- (i-n) 

*J=0,1 x,y£Z d 

Note that (1.1) is the translation invariant convolution equation and admits a simple 
structure in the Fourier space. As a consequence, Fourier representation plays a central role 
in our proofs of properties I and II. On the other hand, Fourier transform alone does not 
suffice in proving III, since our main condition (1.7) is stated in the coordinate space and its 
equivalent interpretation in Fourier space is obscure. 

Property I follows by the method [22]: we prove a uniform bound for the covariance of 
fi t and refer to the Prokhorov Theorem. Property II is deduced from an analysis of the 
oscillatory integral representation of the correlation function in Fourier space. An important 
role is attributed to Lemma 3.1 reflecting the properties of the Fourier transformed correlation 
functions which is derived from the mixing condition. To prove III we exploit the dispersive 
properties of the dynamics (1.1) in coordinate space. The dispersion follows from a stationary 
phase method applied to the oscillatory integral representation of the Green's function in 
Fourier space. The dispersion allows us to represent the solution as a sum of weakly dependent 
random variables by the Bernstein-type 'room-corridor' partition. 

Let us explain in more detail the main idea for the proof of III. First let us consider the 
case n — 1 and the nearest neighbor crystal for which the potential energy has the form 

I £ (y( x -y) u {y)Mxj) = \ E(EK* + e*)-«(:r)| 2 + m 2 Kx)| 2 ), (1.12) 

x,y&Z d xGZ d i=l 
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where m > and e« = (6n, . . . Sid) ■ The solution is represented through the Green's function, 
G(t,x) , 

Y(x,t)= ]T g(t,x-y)Y (y). (1.13) 
yez d 

The long-time asymptotics of the Green's function is analyzed by stationary phase method 
based on the dispersion relation 

d 

u(P) :=V l l 2 {6) = (2]T(l-cos^)+m 2 ) 1/2 , 9 E T d , (1.14) 

3=1 

where T d is the real d -torus and V{9) stands for the Fourier transform of V(x). The main 
features of uj for m > are 

i)u(6)^0, 9eT d , and ii) mesC = 0, (1.15) 

where C is the critical set {9 E T d : det Hess u;(0) = 0} and 'mes' stands for the Lebesgue 
measure in T d . The Green's function has distinct asymptotic behavior in three zones of 
(x, t) -space: inside resp. outside the light cone and in the 'buffer zone', which is a small 
conical neighborhood of the boundary of the light cone. The light cone is determined by the 
group velocities Vu;(#) of the phonons, and its boundary is determined by the group velocities 
Vui{9) with "critical" 9 EC , since they correspond to the maximal values of |Vu;(#)| with 
a fixed direction of Voo(9) (cf. (1.16)). Therefore, the buffer zone is determined by the 
velocities Vuj(9) with the 9 from a small neighborhood of the critical set C . The Green's 
function decays rapidly outside the light cone, as t~ d l 2 inside the light cone except for the 
buffer zone, and more slowly in the buffer zone, cf. (1.18). 

Now let us discuss the general case when n > 1 . For n > 1 an additional important 
feature occurs. In this case we have n dispersion relations oj k {9) , k — 1, n , which are the 
eigenvalues of the matrix V 1 ^ 2 (9) . Thus there can be "crossing points" where two or more 
dispersion relations Uk{9) coincide which implies that they are not differentiable, in general. 
In this case the decay of the Green's function generally is slower than t~ d ^ 2 everywhere in 
(x, t) -space. We estimate the decay by the stationary phase method, hence we need smooth 
branches of the dispersion relations Uk(9) at least locally in 9 . We establish the existence 
of the branches outside a set of the Lebesgue measure zero in T d (see Lemma 2.2). For the 
proof we use the advanced variant of the Weierstrass Preparation Theorem from [15] and the 
analytic stratification of analytic sets [12]. 

For n > 1 we define the critical set C as the subset of T d which is the union over 
k = 1, ...,n of all the points 9 either with a nondifferentiable Uk{9) , or with a degenerate 
Hessian of oo k (9) , or with uo k (9) = . Lemmas 2.2, 2.3 imply that mes C = which plays the 
central role in all proofs in the paper. The critical set is never empty. For example, let us fix 
k = 1, ...,n and consider the point 9 <E T d with the maximal group velocity |Va;fc(0)| > . 
Then det Hess Uk(9) = since Hessu;*:^) Vu)k{9) = 0: 

fu fa\Y7 fa\\ ^ d 2 u k (9) du k (9) 18^*8^(9)2 , , v 

(Hess^V^^g = 2 89^^9- = °' * = W ' (L16) 

3 % J 3 1 j 3 

provided the derivatives exist. Thus even for d — n — 1 the uniform in x G R decay of the 
Green's function is slower than t" 1 ^ 2 since u"(9) vanishes in some points. To overcome this 
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difficulty, in [2] it is required that uj"'{6) 7^ at points with ui"{6) = . Then the uniform 
decay of the Green's function is t -1 / 3 which suffices in the case d — 1 together with an 
additional assumption on the higher moments of the initial measure. In contrast the critical 
set and the slow decay of the Green's function do not occur for the Klein-Gordon equation 
analyzed in [4, 6]. 

For d, n > 1 Suhov and Shuhov have proved in [19] the convergence of the covariance, 
(1.9), for a simple singularity of uj k {6) (in Arnold's terminology [1]) in the points 9 G C 
with the degenerate Hessian. However, similar detailed analysis of all degenerate points for 
d, n > 1 seems to be impossible. We avoid it by a novel 'cutoff' strategy which allows us 
to cover the general case when the Lebesgue measure of the critical set C is zero. Namely, 
we choose an e > and split the Fourier transform of the solution in two components 
Y(9,t) = Yf(9,t) + Y g (9,t) where Yf(9,t) = outside the e -neighborhood of the critical set 
C while Y g (9,t) = inside the e/2 -neighborhood of C . First, we use the mixing condition 
to estimate the contribution from the 'critical' component Yf : we prove that it is small in the 
mean, i.e. its dispersion is negligible uniformly in t > , if e > is sufficiently small. This 
follows from the identity mes C = since the Fourier transforms of the initial correlation 
functions are absolutely continuous due to the mixing condition. A further step is to develop 
a Bernstein type argument to prove the Gaussian limit for the main 'noncritical' component 
Y g . We write it in the form (1.13): 

Y g (x,t)= J2 g (t,x-y)Y (y) (1.17) 
y ez d 

where G g (t, x—y) is the 'truncated' Green's function which is defined similarly to Y g (x, t) : its 
Fourier transform Q g (t, 9) is zero inside the e/2 -neighborhood of C . Then all the dispersion 
relations Uk(0) are smooth and nondegenerate on the support of G g (t, 9) , hence the truncated 
Green's function has the standard decay 



G g (t,x-y) < < 



' Ct~ d / 2 , \y-x\< ct 

(1.18) 

[ C p (\t\ + \x-y\ + iy p , \y-x\> ct 



with some c > and any p > , cf . (5.2), (5.3). Therefore, the representation (1.17) 
demonstrates that for a fixed x G Z d , the main contribution to Y g (x,t) comes from the 
section B t {x) = {y G Z d : \y — x\ < ct} of the light cone at time t . The "volume" of the 
section (i.e. the number of the points y G Z d n B t (x) ) is |S t (x)| ~ t d . Therefore, (1.17) 
becomes roughly speaking, 

E Y (y) 

Y g (x,t) ~ veBt(x) ^ , t^oo. (1.19) 

I -By- 1 



This implies the Gaussian limit by the Ibragimov-Linnik Central Limit Theorem [13], since 
the random values Y (y) are weakly dependent because of the mixing condition (1.7). 

Remarks 1.1 i) Physically, the asymptotics (1.18) reflects the isotropic phonons propaga- 
tion of phonons in the noncritical spectrum. The isotropy provides a 'dynamical mixing' 
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which leads to the Gaussian behavior by the statistical mixing condition (1.7). So the con- 
vergence to the statistical equilibrium (1.8) is provided by the both kinds of the mixing 
simultaneously: the statistical mixing condition (1.7) and the dynamical mixing (1.18). 

ii) The degree —d/2 in (1.18) is related to the energy conservation since the Hamiltonian 
(1.4) is a quadratic form. Roughly speaking, (1.18) means the 'energy diffusion', and the 
degree —d/2 resembles the diffusion kernel. 

Finally, let us comment on our conditions concerning the interaction matrix V(x) . We 
assume the conditions E1-E4 below which in a similar form appear also in [2, 14]. El means 
the exponential space-decay of the interaction in the crystal. E2 resp. E3 means that the 
potential energy is real resp. nonnegative. E4 eliminates the constant part of the spectrum 
and ensures that mesC = (cf. (1.15)). We also introduce a new simple condition E5 for the 
case n > 1 which eliminates the discrete part of the spectrum for the covariance dynamics. 
It can be considerably weakened to the condition E5' from Remark 2.10 Hi) . For example, 
the condition E5' holds for the canonical Gaussian measures which are considered in [14]. 
We show that the conditions E4 and E5 hold for "almost all" matrix- functions V(-) with 
the finite range of the interaction. 

Furthermore, we do not require that uj k (6) ^ 0, 9 G T d : note that lo(0) = for the 
elastic lattice (1.14) in the case m — . Our results hold whenever mes {6 G T d : Uk{0) = 0} = 
. To cover this case we impose the new condition ES which is roughly speaking necessary 
and sufficient for the uniform bounds of the covariance. It can be simplified to the stronger 
condition 

\\V- 1 {e)\\eL 1 (T d ) (1.20) 

from [14], which holds for the elastic lattice (1.14) if either n > 3 or m > 0. The condition 
(1.20) is equivalent to ES for the canonical Gibbs measures considered in [14]. However, 
(1.20) does not hold in some particular interesting cases: for instance, for the elastic lattice 
(1.14) in the case n — 1, 2 and m = , as it is pointed out in [14]. 

The main results of our paper are stated in Section 2: Theorem A in Section 2.4, and its 
application in Section 2.5.4. The convergence (1.9) and the compactness I are established in 
Section 3, and the convergence (1.10) in Sections 4 to 8. Section 9 concerns the ergodicity 
and the mixing properties of the limit measure. In Appendix we analyze the crossing points 
of the dispersion relations. 

2 Main results 
2.1 Dynamics 

We assume that the initial date Y belongs to the phase space H a , a G R , defined below. 

Definition 2.1 H a is the Hilbert space of pairs Y = (u(x),v(x)) of R™ -valued functions 
of x G Z d endowed with the norm 

\\Y\\l = E (K*)| 2 + Hx)\ 2 )(l + M 2 )* < oo. (2.1) 

x£Z d 
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We impose the following conditions El - E5 on the matrix V . 

El There exist constants C, a > such that |Vm(z)| < Ce^ a ^, k,l G I„ := {1, ■■■,n}, z G 
Z d . 

Let us denote by V{9) := (v k i(9)) u , y , where V kl {9) = E V w (z)e <a * 9 G T d , and T rf 

\ / fe, (El n zgZ d 

denotes the d -torus T d = K d /2nZ d . 

E2 V is real and symmetric, i.e. Vik(-z) = Vki(z) G R, k, I G I n , -2 G Z d . 

The condition implies that is a real-analytic Hermitian matrix-function in 9 G T d . 

E3 The matrix V(9) is non-negative definite for each 9 G T d . 

The condition means that the Eqn (1.1) is hyperbolic like wave and Klein-Gordon Eqns 
considered in [6] - [8]. Let us define the Hermitian non- negative definite matrix 

£1(6) := (V(9)) 1/2 > (2.2) 

with the eigenvalues Uk(9) > 0, k G I n , the dispersion relations. For each 9 G T d the 
Hermitian matrix Q(0) has the diagonal form in the basis of the orthogonal eigenvectors 
{e k (9):kel n }: 

1^(9) ... \ 
Q(9) = B(9) o ••. B*(9), (2.3) 
V ... co n (9) j 

where B(9) is a unitary matrix. It is well known that the functions ujk(9) and B(9) are 
real-analytic outside the set of the 'crossing' points 9* : ^(0*) = uoi(9*) for some I ^ k. 
However, generally the functions are not smooth at the crossing points if 0Jk(9) ^ oji(9) . 
Therefore, we need the following lemma which we prove in Appendix (cf. [21, Lemma 1.1 ]). 

Lemma 2.2 Let the conditions El, E2 hold. Then there exists a closed subset C T d s.t. 

i) the Lebesgue measure of is zero: 

mesC* = 0. (2.4) 

ii) For any point G T d \C* there exists a neighborhood 0{Q) s.t. each dispersion relation 
u)k{9) and the matrix B{9) can be chosen as the real-analytic functions in 0(Q) . 

Hi) The eigenvalues uik{9) have constant multiplicity in T d \ C* , i.e. it is possible to enu- 
merate them so that we have for 9 G T d \ C* , 

^(9) = . . .=w n (6), uj ri+1 (9)= ... =LJ r2 (9), . . . , v rB+ i(6) = . . .=uj n (9), (2.5) 

u ra (0) ^ u ru (9) if ff^i/, 1 < TV, r u < r s+1 := n. (2.6) 
iv) The spectral decomposition holds, 

£1(9)^^(9)^(9), 0eT d \C*, (2.7) 
i 

where H a (0) is the orthogonal projection in R n which is real- analytic function of 9 G T d \C* . 
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Below we denote by u k (9) the local real-analytic functions from Lemma 2.2 ii) . Our next 
condition is the following: 

E4 D k {9) ^ , Vfc G I n , where D k {9) := det (f^®)* , # G T d \ . 

Let us denote C := G T d : det \>(0) = 0} and C k := {9 G T d \ : D k (9) = 0} , fc G I n . 
The following lemma is also proved in Appendix. 

Lemma 2.3 Let the conditions El - E4 hold. Then mes C k = for k = 0, 1, ...,n . 
Our last condition on V is the following: 

E5 For each k ^ I the identity w k (9) — ivi(0) = const_ , G T d does not hold with 
const _ 7^ , and the identity c^fc(^) + tvi(9) = const + does not hold with const + ^ . 

This condition holds trivially in the case n = 1 . 

We show that the conditions E4 and E5 hold for "almost all" functions V satisfying 
the conditions El, E2. More precisely, let us fix an arbitrary N > 1 and denote by TZn 
the set of the "finite range" interaction matrices V with V(x) = for maxj \xi\ > N , and 
satisfying the condition E2. In Appendix we prove the following lemma. 

Lemma 2.4 For any N > 1 the conditions E4 and E5 hold for the matrix-functions V 
from an open and dense subset of 1Zn ■ 

The following proposition is proved in [14, p. 150] and [2, p. 128]. 

Proposition 2.5 Let El and E2 hold, and a G R . Then 

i) for any Y G H a there exists a unique solution Y{t) G C(R, 7i a ) to the Cauchy problem 
(1.2). 

ii) The operator U(t) : Y \— > Y{t) is continuous in H a . 
Proof Applying the Fourier transform to (1.2), we obtain 

Y(9,t) = A(9)Y(9,t), t G R, Y(0) = Y , (2.8) 

where 

A ^-( J m I)- 9eTi < 2 - 9 > 

Note that Y(-,t) G D'(T d ) for t G R. On the other hand, V"(0) is a smooth function by 
El. Therefore, the solution Y(9,t) of (2.8) exists, is unique and admits the representation 
Y{9,i) = exp (^A(9)tjY (9) . It becomes (1.13) in the coordinate space, where the Green's 
function Q (t, z) admits the Fourier representation 

G(t, z) := F e -i 2 [exp (A(9)t)} = (2n)- d J e^ 6 exp (A(9)t)d9. (2.10) 

Hence, by the partial integration, Q(t, z) ~ \z\~ p as \z\ — > oo for any p > and bounded 
\t\ since ^4(0) is the smooth function of 9 G T d . Therefore, the convolution representation 
(1.13) implies Y(t) G H a . □ 
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2.2 The convergence to statistical equilibrium 

Let (Q, X, P) be a probability space with expectation E and let B(7i a ) denote the Borel 
a -algebra in H a . We assume that Y = Y (cu, •) in (1.2) is a measurable random function 
with values in (H a , B(H a )) . In other words, for each x G Z d the map uo i— > Y^(a;,a;) is 
a measurable map f2 — > R 2n with respect to the (completed) a -algebras X and £>(R 2n ) . 
Then F(t) = U(t)Y is again a measurable random function with values in (Ha,B(H a )) 
owing to Proposition 2.5. We denote by fio(dYo) a Borel probability measure on 7i a giving 
the distribution of the Y . Without loss of generality we assume (f2, X, P) = (H a , B(H a ), /io) 
and Y (u,x) = uj(x) for fio(du) -almost all u G H a and each x G Z d . 

Definition 2.6 fi t is a Borel probability measure in H a which gives the distribution of 
Y(t) : 

m{B) = no{U(-t)B), MB G B(H a ), teR. (2.11) 

Our main goal is to derive the convergence of the measures ji t as t — > 00 . We establish 
the weak convergence of fi t in the Hilbert spaces H a with a < —d/2 : 

l^t — " Moo as t —> 00, (2.12) 
where jj,^ is a limit measure on the space TL a , a < —d/2 . This means the convergence 

J f(Y)nt(dY) - I /(rVo^n * - °o, (2.13) 
for any bounded continuous functional / on 7i a . 

Definition 2.7 TTie correlation functions of the measure ji t are defined by 

Qi j (x,y) = E(Y i (x,t)®Yi(y,t)), i,j = 0,1, x,yeZ d , (2.14) 

if the expectations in the RHS are finite. Here Y l (x,t) are the components of the random 
solution Y(t) = (F (-,t),F 1 (-,t)) . 

For a probability measure /i on H a we denote by ft the characteristic functional (Fourier 
transform) 

/}(*) = J exp(i(Y, *)) /i(rfF), $GD. 
A measure /i is called Gaussian (of zero mean) if its characteristic functional has the form 

/i(*) = exp{-^Q(*,*)}, (2.15) 

where Q is a real nonnegative quadratic form in T> . A measure fi is called translation- 
invariant if fi(T h B) = /2(B), B G B(H a ), h G Z d , where T h Y(x) = Y(x-h), x G Z d . 
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2.3 The mixing condition 

Let 0(r) denote the set of all pairs of the subsets A, B C Z d at distance dist (A, B) > r and 
cr(A) be the a -algebra in 7i a generated by Y(x) with x G A. Define the Ibragimov-Linnik 
mixing coefficient of a probability measure /i on Tt a by (cf. [13, Definition 17.2.2]) 

\fi (AnB)-fi (A)fi (B)\ 
<p(r) := sup sup — . (2.16) 

(AB)eO(r) A G (7(^,5 G (7(B) 

/io(5) > 

Definition 2.8 The measure /i satisfies strong, uniform Ibragimov-Linnik mixing condition 
if <f(r) ^ as r — > oo. 

Below, we specify the rate of decay of (p (see condition S3). 



2.4 Statistical conditions and results 

We assume that the initial measure /x satisfies the following conditions S0-S3: 

50 /io has zero expectation value, EY (x) = , xeZf 

51 /io has translation-invariant correlation matrices, i.e. Eqn (1.5) holds for x, y G Z d . 

52 fj, has a finite mean energy density, i.e. Eqn (1.6) holds. 

53 /io satisfies the strong uniform Ibragimov-Linnik mixing condition with 

Tp:= J r d ~ V /2 (r) dr < oo. (2.17) 
o 

We will deduce from SO - S3 that gj j G C(T d ) , i, j = 0, 1 (see Lemma 3.1). This makes 
sense of our last condition ES concerning the initial covariance and the matrix Q(6) . We 
need it only in the case when Co ^ i.e. det V(9) — for some points 9 G T d : 

ES 11^(0)^(0)^(0)11 G L\T d ) for i,j = 0,1. 

This condition follows from SO - S3 if i — j — or C = . 



Next introduce the correlation matrix of the limit measure /ioo • It is translation-invariant 
(cf. (1.5)): 

Q 00 (x,y)^(qg(x-y)). j=QX (2.18) 

In the Fourier transform we have locally outside the critical set C* (see Lemma 2.2) 

q£(9) = B(9)M£(9)B*(9), i,j = 0,1, (2.19) 

where B(9) is the smooth unitary matrix from Lemma 2.2 ii) and M^(9) is n x n -matrix 
with the smooth entries (m£(0)) = Xki(B*(9)M l Q j (9)B(9)) . Here we set (see (2.5)) 

_/l if fc,Z G (tv_i,7v], o- = l,...,s + 1, , v 

XM " \ otherwise (2 ' 20) 
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with r := , r s+ i := n , and 

/ q° °(9) + n-\9) ql\6) Q~\9) <%\9) - n~\9) q™{9) Q(9) \ 

The local representation (2.19) can be expressed globally as 



M (9) := \ 



(2.21) 



q£(9) = S f2n*(9)Mi ) j (9)n a (9), 9eT d \C*, i,j = 0,1, (2.22) 

a=l 

where H a (9) is the spectral projection from (2.7). 

Remark 2.9 The condition ES implies that (M l j ) kl G L\T d ) , k, I G I„ . Therefore, (2.22) 
and (2.4) imply that also G L\T d ) , k,l G I n . 

Theorem A Let d, n > 1 , a < —d/2 and assume that the conditions El— E5 ; SO— S3 hold. 
If C 7^ , then we assume also that ES holds. Then 

i) the convergence in (2.12) holds. 

ii) The limit measure ji^ is a Gaussian translation-invariant measure on TL a . 
Hi) The characteristic functional of is the Gaussian 

=exp{-^Q 0O (*,*)}, *GD, (2.23) 

where is the quadratic form defined in (1.11). 

iv) The measure /z^ is invariant, i.e. [U(t)]*/j,oo = Hoo , t G R. 

Remarks 2.10 i) In the case n — 1 , the formulas (2.21), (2.22) become 

/ «oo I . ,-2 -11 -01 -10 \ 

I % + UJ % % - % 



-10 -01 -11 , , ,2 -00 / 



ii) The uniform Rosenblatt mixing condition [18] also suffices, together with a higher 
power > 2 in the bound (1.6): there exists 5 > such that 

e(\u (x)\ 2+s + \v (x)\ 2+s ) <oo. 

r+oo 

Then (2.17) requires a modification: / r d ~ 1 a p (r)dr < oo , with p = mii^^r^, |) , where 

a(r) is the Rosenblatt mixing coefficient defined as in (2.16) but without Ho(B) in the 
denominator. With these modifications, the statements of Theorem A and their proofs 
remain essentially unchanged. 

Hi) The arguments with condition E5 in Proposition 3.2 (see (3.7)-(3.13) below) demon- 
strate that the condition could be considerably weakened. Namely, it suffices to assume 

E5' If for some k ^ I we have either uj k (9)+uji(9) = const + ^ or oo k (9)—uji(9) = const_ ^ 
, then 

(B\9)q i j{9)B{9)) ki = 0, 9 G T d , i,j = 0, 1. (2.24) 
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The assertions i) — in) of Theorem A follow from Propositions 2.11 and 2.12: 

Proposition 2.11 The family of the measures {/i t , t G R} is weakly compact in H a with 
any a < —d/2, and the bounds hold: 

sup£||l/(*)yo||a<°o- (2-25) 

t>o 

Proposition 2.12 For every \& G T> , the convergence (1.10) holds. 

Proposition 2.11 ensures the existence of the limit measures of the family {fi t , t G R} , while 
Proposition 2.12 provides the uniqueness. Propositions 2.11 and 2.12 are proved in Sections 
3 and 4-8, respectively. 

Theorem A iv) follows from (2.12) since the group U(t) is continuous in H a by Propo- 
sition 2.5 ii) . 

2.5 Examples and applications 

Let us give the examples of the equations (1.1) and measures /x which satisfy all conditions 
E1-E5, S0-S3 and ES. 

2.5.1 Harmonic crystals 

All conditions E1-E5 hold for ID crystal with n = 1 considered in [2]. For any d > 1 and 
n — 1 consider the simple elastic lattice corresponding to the quadratic form (1.12) with 
m ^ . Then V(x) = F^ x uj 2 (6) with u(9) defined by (1.14), satisfies E1-E4 with C* = . 
In these examples the set C is empty, hence the condition ES is superfluous. Condition E5 
holds trivially since n = 1 . 

2.5.2 Gaussian measures 

We consider n = 1 and construct Gaussian initial measures /i satisfying S0-S3. We will 
define /i by the correlation functions (£o{x — y) which are zero for i ^ j , while for i — 0, 1 , 

$(6) := F z _ e q^z) G L\T d ), ${6) > 0. (2.26) 

Then by the Minlos theorem [11] there exists a unique Borel Gaussian measure /io on H a , 
a < —d/2, with the correlation functions qo(x — y) . The measure /xo satisfies S0-S2. 
Further, let us provide, in addition to (2.26), that 

0, \z\>r . (2.27) 

Then the mixing condition S3 follows with ip(r) = , r > r , since for Gaussian random 
values the orthogonality implies the independence. For example, (2.26) and (2.27) hold if we 
set q l l (z) = f(z 1 )f(z 2 ) ■ ■■ ■ ■ f(z d ) , where f(z) = v - \z\ for \z\ < u and f(z) = for 
\z\ > uq with uq := [r /^/d] (the integer part). Then by the direct calculation we obtain 

f(9) = 1 _ cosQ , e T 1 , and (2.26) holds. The measure /io is nontrivial if ro > Vd: 

otherwise u = , so q % o(z) = 0, and the measure /i (dY ) is concentrated at the point 
Y = 0. 



12 



2.5.3 Non-Gaussian measures 



Let us choose some odd bounded nonconstant functions f°, f 1 G C(R) and consider a ran- 
dom function (Y°(x), Y l (x)) with the Gaussian distribution /x from the previous example. 
Let us define /x|$ as the distribution of the random function (f°(Y°(x)), f 1 (Y 1 (x))) . Then 
SO - S3 hold for /x|$ with corresponding mixing coefficients ip*(r) = for r > r . The 
measure /Xg is not Gaussian if the functions f° , f 1 are bounded and nonconstant. 

2.5.4 From statistical chaos to the Gibbs measure 

Let us consider the initial measures which satisfy SO - S3, and with the correlation functions 

(q$) a (x-y):=E(YZ(x,0)Y?(y,0))=T i 6 ij 6 kl 6 xy , i,j = 0,1, k, I G I„, x,x/GZ d , (2.28) 

where T 0) i > . These correlations correspond to the "chaos" with the zero correlation 
radius and noncorrelated components. Such measures exist on H a with a < —d/2 by the 
Minlos Theorem [11]: for example, the "white noise" which is the corresponding Gaussian 
measure. Let us consider the crystal satisfying the conditions El - E4 and (1.20). Then also 
the conditions E5', ES hold, so Theorem A is applicable (see Remark 2.10 Hi) ) : it implies 
the convergence (2.12) to the Gaussian measure /Xqo with the covariance (2.21), (2.22). 

Additionally, let us assume that T = which physically means that only the initial 
velocities contribute, and initial deviations are adjusted to zero. Then the formulas (2.21), 
(2.22) become 

/ V~\e) \ 



q oo (6) = M (6) = ^ 







(2.29) 



V u V kl Jk,iei n / 



According to (1.3), this means that the limit measure ji^ coincides with the Gibbs canonical 
measure corresponding to the temperature ~ T\ . In a more general framework, the limit 
measure is close to the Gibbs measure if the radius of the initial correlations is small in a 
suitable scaling limit (cf. [6, Proposition 4.2]). 

3 Convergence of covariance and compactness 
3.1 Mixing condition in terms of spectral density 

The next Lemma reflects the mixing property in the Fourier transforms of initial corre- 
lation functions $ . Condition S2 implies that qo(z) is a bounded function. Therefore, its 
Fourier transform generally belongs to the Schwartz space of tempered distributions. 

Lemma 3.1 Let the conditions SO - S3 hold. Then $ e C(T d ) , i, j = 0, 1 . 

Proof It suffices to prove that 

tj${z) e l\z d ). (3.1) 

Conditions SO - S3 imply by [13, Lemma 17.2.3] (or Lemma 8.2 i) below): 

|go 5 WI<cW /2 (M), zez d , (3.2) 
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where eo is defined by (1.6). Therefore, (2.17) implies (3.1): 

E |g^WI<CeoE^ 1/2 (kl)<oo. □ 

z£Z d zeZ d 

3.2 Oscillatory integral arguments 

In this section we uniformly estimate and check the convergence of the correlation matrices of 
measures /i t with the help of the Fourier transform. The condition SI and the translation- 
invariant dynamics (1.1) imply that 

Q?(x,y) = J Y\x) ® Yi{y)m{dY) = ${x - y), x,y G Z d . (3.3) 

Proposition 3.2 i) The correlation matrices q t 3 {z), i,j = 0,1, are uniformly bounded 

sup sup |?t" ? (^)| < oo. (3.4) 
t>o zez d 

ii) The correlation matrices ql J (z) , i,j — 0, 1 , converge for each z e Z d , and 

q?(z)^q£(z), t^oo, (3.5) 
where the functions q 1 ^ (z) are defined above. 

Proof. For brevity, we prove (3.4) and (3.5) for % = j = . In all other cases the proof of 
(3.5) is similar. The solution to the Cauchy problem (1.1) is 

u(x,t) = (2n)- d J e- lx - e (costtt Y°(6) + smttt IT 1 ^))^, 

where Q = Q(6) is the non-negative definite Hermitian matrix defined by (2.2). Furthermore, 
the translation invariance (1.5) implies that 

e(y*{0) ® Yi{9')) = (2n) d 8(e + tf)${9), i,j = 0, 1. (3.6) 

Hence, 

q™{x-y) := E (u(x , t) ® u(y , t)) 

= {2n)- d j e ~ l9 ^ [cos Qt q° °(6) cos fit + sin fit cos fit 

+ cos fit gg 1 ^" 1 sin n* + sin fi* fi" 1 ^ 1 ^)^ 1 sin Qt]d9. (3.7) 
Therefore, the bound (3.4) with i — j — follows from Lemma 3.1 or condition ES if 

Let us check that the convergence (3.5) with i — j — also follows since the oscillatory 
integrals in (3.7) tend to zero. Consider for example the last term in the integrand of (3.7). 
We rewrite it using (2.3), in the form 

Lg(6,t) := sinfit ST 1 ft 1 (6) fi _1 sinfit = B(0)(smu k t A%(6) sinc^t) B*{6), (3.8) 
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where A u (9) := B*(9)Q~ 1 ql 1 (9)Q~ 1 B(9) . However, at this moment we have to choose 
certain smooth branches of the functions B{9) and Uk{9) since we are going to apply the 
stationary phase arguments which require a smoothness in 9 . To make it correctly, we cut 
off all singularities. First, we define the combined critical set 

C := UfcCfeUC* UC . (3.9) 

Then Lemmas 2.2, 2.3 imply the following lemma. 

Lemma 3.3 Let conditions El - E4 hold. Then mes C = . 

Second, fix an e > and choose a finite partition of unity 

M 

f(9) + g(9) = 1, g{9) = £ g m (6), 9 e T d , (3.10) 

m=l 

where /, g m are nonnegative functions from C£°(T d ) , the supports of g m are sufficiently 
small and 

supp/ C {9 e T d : dist(0,C) < e}, supp g m C {9 e T d : dist(0,C) > e/2}. (3.11) 
Now (3.8) can be rewritten as 

Ll\e,t) = f(e)Ll\e,t) 



+\Y,3rn{9)B{9)Ucos{u k -u l )t - cosicuk + cu^AlKe)) B*(9). (3.12) 

By Lemma 2.2 and the compactness arguments, we can choose the supports of g m so small 
that the eigenvalues ujk(9) and the matrix B(9) are real-analytic functions inside the suppg m 
for every m : we do not mark the functions by the index m to not overburden the notations. 

Let us substitute (3.12) into the last term of (3.7) and analyze the Fourier integrals with 
/ and g m separately. The integral with / converges to zero uniformly in t > , as e — > . 
Indeed, by Lemma 3.3 we have 



f(9)L 1 1 (9,t)d9\<C J \\n- 1 (9)q 1 1 (9)n- 1 (9)\\d9 ^0, e^O 



j, '" ' now ■ 1 

T d dist(e,C)<£ 



since the integrand is summable by Lemma 3.1 or condition ES if Co ^ • 

Below we will prove the convergence for the integrals with g m . We will deduce the 
convergence from the fact that the identities oo k (9) ±ui(9) = const ± with the const ± ^ are 
impossible by the condition E5. Furthermore, the oscillatory integrals with ujj,±uji(9) ^ const 
vanish as t — > oo . Hence, only the integrals with U0k{9) — u)i(9) = contribute to the limit 
since uik{9) + uii{9) = would imply 0Jk{9) = 0Ji{9) = which is impossible by E4. Similar 
analysis of the three remaining terms in the integrand of (3.7) gives, 



q™{x-y) = (2n)- d fe- ie ^f(9)Ll 1 (9,t)d9 
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d9 



= (2rr)- d J e~ w{x - y) f (6)1^(6, t) d9 + (2vr)- d | g(9)e- ie ^q^(9) d9 + . . . (3.13) 

according to the notations (2.18) - (2.21). Here A 00 (9) := B*(9)$°(9)B(9) and ' ... ' stands 
for the oscillatory integrals which contain cos(a;fc(#) ±u>i(9))t and sin(u>k(9) ±u>i(9))t with 
u>k{9) ±u>i(9) ^ const. 

The oscillatory integrals converge to zero by the Lebesgue-Riemann Theorem since all 
the integrands in ' ... ' are summable and V(uJk{9) ± uji(9)) = only on the set of the 
Lebesgue measure zero. The summability follows from Lemma 3.1 or the condition ES 
since the matrices B*(9) are unitary. The zero measure follows similarly to (2.4) since 
uik{9) ±u>i(9) ^ const. 

At last, let us prove the convergence (3.5) with % — j — . From the last line of (3.13) 
we know that qf°(x — y) is close to the integral with g if e > is small and t > is 
large. Therefore, the limit of q^°(x — y) as t — > oo coincides with the limit of the integral as 
e — > . Finally, this limit coincides with q^(x — y) since G L 1 (T d ) by Remark 2.9. □ 

3.3 Compactness of measures family 

Proof of Proposition 2.11 The compactness of the measures family {fi t , t G R} will follow 
from the bounds (2.25) by the Prokhorov Theorem [22, Lemma II. 3.1] using the method of 
[22, Theorem XII. 5. 2] since the embedding H a C Hp is compact if a > (3. 

First, the translation invariance (3.3) and Proposition 3.2 i) imply that for x G Z d we 
have 

e t := J [\u (x)\ 2 + \v (x)\ 2 } n t (dY ) = trg, 00 (0) + tr % n (0) < e < oo, * > 0. (3.14) 

Hence by the definition (2.1) we get for any a < —d/2 : 

E\\U(t)Y \\ 2 a = et J2(l + \x\ 2 ) a = C(a)e t <C(a)e<oo, t > 0. □ 
xez d 

4 Duality argument 

To prove Theorem A, it remains to check Proposition 2.12. Let us rewrite (1.10) as follows, 

£expO(y(0,*)}^/}ooW, t^oo. (4.1) 

We will prove it in Sections 5-9. In this section we evaluate (Y(t), by using the following 
duality arguments. Remember that Y G 7i a with a < —d/2 . For t G R introduce a 
'formal adjoint' operator U'(t) from space T> to 7i- a : 

(Y,U'(t)*) = (U(t)Y,V), $GD, YeH a . (4.2) 
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Let us denote by $(•,£) = U'(t)^/ . Then (4.2) can be rewritten as 

(Y(t),^) = (Y ,^(-,t)), tGR. (4.3) 

The adjoint group U'(t) admits the following convenient description. Lemma 4.1 below 
displays that the action of group U'(t) coincides with the action of U(t) , up to the order of 
the components. 

Lemma 4.1 For * = (*°, ^Z 1 ) eV we have 

ty-,t):=U'(t)* = (i,(-,t)M-M ( 4 - 4 ) 
where i/)(x,t) is the solution of Eqn (1.1) with the initial data (u ,v ) = . 

Proof Differentiating (4.2) in t with Y, * E V , we obtain that (Y,U'(t)V) = (U(t)Y,V). 
The group U(t) has the generator A from (1.3). The generator of U'(t) is the conjugate 
operator to A : 

^=(;- v ). (4.5, 

Hence, the representation (4.4) holds with tp(x,t) = — J2 V(x — y)ip(y,t) . □ 

y ez d 

The lemma allows us to construct the oscillatory integral representation for $(x, t) . 
Namely, (4.4), (4.5) imply that in the Fourier representation for $(-,£) = U'(t)^> we have 

4(0, t) = A*(6)$(6, t), $(0, t) = g*(t, 6)4? (6). 
Here we denote (see (2.9)) 

with ft = £1(6) = ft*(0) . Therefore, 

$(x,t) = {2nr d J e- l6x g*(t, 6)4/(6) d6, x e Z d . (4.7) 

Since f(6) + g(6) = 1 by (3.10), we can split $ in two components: 

t) = (2ir)- d J e- t9x g*(t, 6)f(6)V(6) dd + (2n)- d J e~ i9x g*(t, 0)#(0)£(0) dd 

I'd T d 

= $f(x,t) + $ g (x,t), xeZ d , (4.8) 

where each function $/(x, t) and t) admits the representation of type (4.4). By (3.11), 

the Fourier spectrum of $/ is concentrated near the critical set C , while the spectrum of 
$ 3 is separated from C . 
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5 Standard decay in the noncritical spectrum 



We prove the decay of type (1.18) for the 'noncritical' component $ 9 . The function <£> 3 can 
be expanded similarly to (3.12), in the form 

*g(*,t) = E E / 9m{e)e-^^aUeme) 60. (5.1) 
m ±, kei n T d 

By Lemma 2.2 and the compactness arguments, we can choose the eigenvalues ui^iO) and 
the matrices a^(9) as real- analytic functions inside the supp g m for every m: we do not 
mark the functions by the index m to not overburden the notations. 

Lemma 4.1 means that each component 3> g (x, t) , % — 0, 1 , is a solution to Eqn (1.1). To 
prove (4.1), we analyze the radiative properties of $ g (x, i) in all directions. For this purpose, 
we apply the stationary phase method to the oscillatory integral (5.1) along the rays x = v t , 
t > . Then the phase becomes {9v ± Uk(9))t , and its stationary points are the solutions to 
the equations v = ^fVuk(9) . We collect all necessary asymptotics in the following lemma 
(cf. (1.18)). 

Lemma 5.1 For any fixed ^ G V and g(9) G C£°(T d \ C) the following bounds hold: 

i) sup \<$> g (x,t)\ < C r d/2 . (5.2) 

x£Z d 

ii) For any p > there exist C p ,^ g > s.t. 

\$ g (x,t)\ <C p (\t\ + \x\ + l)- p , \x\> lg t. (5.3) 

Proof Consider $ s (x, t) along each ray x = vt with arbitrary v e R d . Substituting to 
(5.1), we get 

* 9 (vt,t) = E E / fc^^^) d9. (5.4) 

This is a sum of oscillatory integrals with the phase functions 0^ (9) — 9v ± c^fc(^) an d 
the amplitudes a^(^) which are real-analytic functions of the 9 inside the suppg m . Since 
u>k(9) is real- analytic, each function 0^ has no more than a finite number of stationary 
points 9 G suppg m , solutions to the equation v = =pVuk(9) . The stationary points are 
non-degenerate for 9 G supp g m by (3.11), (3.9) and E4 since 

det (^^) = ±Dk ^ ^ °' G SU PP^- ( 5 - 5 ) 

At last, ^(9) is smooth since f e D. Therefore, § g (vt,t) = 0{tr d l 2 ) according to the 
standard stationary phase method [10, 17]. This implies the bounds (5.2) in each cone 
\x\ < ct with any finite c. 

Further, denote by v g := max m max fceIn max |Vc<j fc (#)|. Then for \v\ > v g the sta- 

6esuppg m 

tionary points do not exist on the supp g . Hence, the integration by parts as in [17] yields 
& g (vt,t) = G(t~ p ) for any p > . On the other hand, the integration by parts in (5.1) 
implies similar bound $ g (x,t) = 0((t/\x\) 1 ^ for any I > 0. Therefore, (5.3) follows with 
any 7 9 > v g . Now the bounds (5.2) follow everywhere. □ 
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6 Contribution of critical set 

We are going to prove (4.1). Rewrite it using (4.3): 

Eexp{«(y , ->0, t^oo. (6.1) 

The splitting (4.8) gives (F , $(•, t)) = (y , $/(-,£)) + (F , $ s (-,t)) . Our main argument is 
that the contribution of (Y , $/(•, £)) to (6.1) has a small dispersion. We will deduce this from 
Lemmas 3.1, 3.3. At first, let us estimate the difference in (6.1) by the triangle inequality: 

|£exp{i(y ,$M))}-£oo(*)l 
< |^exp{i(^)yo,*>}-Sexp{i(y ,^(^0>}| + |/ioo(^-/iooWI 

+\E exp{t(Y , *))} - Aoo(^ 9 )| =/ + // + (6.2) 

where \l/ s := -F -1 ^^)^^)] = $ 9 (-,0) . Let us consider each of the three terms separately. 

I. The first term / = I(e,t) represents the contribution of the neighborhood of the critical 
set {9 G T d : dist(6>, C) < e} and tends to zero as e — > uniformly in t > 0. Namely, by 
the Cauchy-Schwartz inequality, 

/ = |£ e <<ro,*(^)> _£ e i<ro,* 9 (-.t)>| < £ | e l < y °, */(■.*)> - 1| < c(£|(Y , $/(-,t))| 2 ) 1/2 . (6.3) 

Using the Parseval identity and (4.8), we get 

E|(y ,$ / (-,t))| 2 = (27r)- 2 ^|(y w,/(^<|.^,t))| 2 

= (27r)- 2d {E(Y (9) ®y^)), f(9)f(9')g*(t,9)^(9) ® g*(t,9')^(9')). (6.4) 

Now take into account that £(y o (0) <g) Y o (0')) = (^) d S(0 - 9')q (9) similarly to (3.6). Then 
(6.4), (4.6), (3.11) and the bounds < f(9) < 1 imply 

E\<y ,* f (-,t))\ 2 <Ci E / \\n- i (9)qi j (9)n-i(9)\\d9^o, e^O 

i ^= ' 1 dist(e,c)< £ 

owing to Lemma 3.3 since the integrand is summable. The summability follows from Lemma 
3.1 or condition ES if C ^ . 

II. The second term II = 11(e) tends to zero as e — > . Indeed, 

Qoo(^,^) = (2tt)- m E J(^(e),g(6)^(e)<»g(e)¥^dB^Q 0O (%^) 1 e^O 

i,j=0 Td 

by the Lebesgue Dominated Convergence Theorem since < g(9) < 1 and q 1 ^ e L 1 (T a! ) by 
Remark 2.9. Hence for the Gaussian measure , we get by (2.23) 

|Aoo(* fl ) - Aoo(*)| = I exp{-^Q 00 (^ 9 , ^ 9 )} - exp{-iQ 00 (^, *)}| -> 0, £ -> 0. 
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III. To prove Proposition 2.12, it remains to check that for any fixed e > , we have 

III(e,t) = \Eexp{i(Y ,%(-,t))}-ji oo (y g )\^0, t -> oo. (6.5) 
We prove (6.5) in Section 8 using the Bernstein arguments of the next section. 

7 Bernstein's 'rooms-corridors' partition 

Our proof of (6.5) is similar to the case of the continuous Klein-Gordon equation in R d [6]: 
all the integrals over H d become the series over Z d etc. Another novelty in the proofs is the 
following: in the case of the Klein-Gordon equation we have &(x, t) = for \x\ >t+ c(^/) , 
while for the discrete crystal we have (5.3) instead. 

Let us introduce a 'room-corridor' partition of the ball {x G Z d : \x\ < 7 9 t} with 7 9 
from (5.3). For t > we choose below A u p t G N (we will specify the asymptotical relations 
between t , A t and p t ) . Let us set h t — A t + p t and 

a? = jht, V = a 3 + A t , j G Z, N t = [{i g t)/h t \. (7.1) 

We call the slabs R{ = {x e Z d : \x\ < N t h t , a j < x d < b>} the 'rooms', C\ = {x G Z d : 
\x\ < N t h t , W < x d < a j+1 } the 'corridors' and L t = {x G Z d : \x\ > N t h t } the 'tails'. Here 
x = (xi, . . . , Xd) , A t is the width of a room, and p t is that of a corridor. Let us denote by 
Xt the indicator of the room R\ , $ that of the corridor C\ , and r\ t that of the tail L t . 
Then 

EM*) + ^W] + ifc(x) = 1, x G Z d , (7.2) 
w t -i 

where the sum J2t stands for ■ Hence we get the following Bernstein's type represen- 

j=-N t 

tat ion: 

fa, *,(-,*)> = Ej< y o,^(-,*)> + (Yo,e t %(;t))] + (lo,^ ff (-,i)». (7.3) 
Let us introduce the random variables r{ , c J t , l t by 

W = (y ,^(-,*)>, ^ = (F ,^(-,t)), z t = (y ,^(-,0>- (7-4) 

Then (7.3) becomes 

O r o,^(-,*)> = E t (W' + ^) + if (7-5) 
Lemma 7.1 Let SO— S3 hold. The following bounds hold for t > 1 : 

£|ri'| 2 < C7(* 9 ) A t /t, Vj, (7.6) 
E\c{\ 2 < C(^ g ) Pt /t, Vj, (7.7) 
E\l t \ 2 < C p (V g ) (l + t)~ p , \/p>0. (7.8) 



Proof We discuss (7.6), and (7.7), (7.8) can be done in a similar way (the proof of (7.8) 
additionally uses (5.3)). Express E\r{ | 2 in the correlation matrices. Definition (7.4) implies 
that 

E\r\\ 2 = (xi(x)xi(y)q (x - y),$ g (x,t) ® $ g {y,t)). (7.9) 
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According to (5.2), Eqn (7.9) implies that 

E\H? < Ct- d ^xi(x)\\qo(x-y)\\ 



x,y 



= cr d £ xi(x) <CA t /t, (7.10) 

X z 

where ||?o(^)|| stands for the norm of a matrix (qo(zf) . Therefore, (7.10) follows as 

IMOIIe'W b y (3-i)- ° 
8 Ibragimov-Linnik Central Limit Theorem 

In this section we prove the convergence (6.5). As was said, we use a version of the Central 
Limit Theorem developed by Ibragimov and Linnik [13]. If Qoo(^ g , = , (6.5) is obvious. 

Indeed, \Ee^'^ -1| < E\ (Y , %(-, t))\ < (e(Y , $,(-, t)> 2 ) 1/2 = (Q t (V g , ^g))^ , where 
Qt(fyg, tyg) — > Qoo(^ g , ^ g) = 0, as t — > oo. Thus, we may assume that for a given $gP, 

Goo(*„* fl )^0. (8.1) 

Let us choose < 5 < 1 and 

Pt-t 1 - 5 , A t ~-^, t^oo. (8.2) 
logt 

Lemma 8.1 The following limit holds true: 

^(^) + (f) 1/2 )+^ 2 (^ 1/2 (p*) + f)-0, f-oo. (8.3) 
Proof The function <p(r) is non-increasing, hence by (2.17), 

r r 

rV /2 (r) = d y" s d -V 1/2 (0 rfs < rf| s d -y /2 (s) < Cy? < oo. (8.4) 



Furthermore, (8.2) implies that h t = A t + p t ~ , t — > oo . Therefore, iV* ~ — ~ logt . 

logt ^ 

Then (8.3) follows by (8.4) and (8.2). □ 
Proof of (6.5) By the triangle inequality, 

\Eew{i(Y ,* g (-,t))} < \E exp{t(Y ,$ g (-,t))} - E exp{tJ2 t rl}\ 

+ 1 expi-^i^'l 2 } - exp{-iQ 00 (^, *,)}| 

+\Eexp{i%2 t ri} - exp{-^ t ^IW | 2 }| 
= I 1 + I 2 + I 3 . (8.5) 
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We are going to show that all the summands l\ , h , ^3 tend to zero as t — > 00 . 
S'tep z) Eqn (7.5) implies 

h = \EeM^2/t}(eM^2 t 4 + ik}-l)\ 

< Cj2 t E\4\ + E\l t \ < C^ t (E\4\ 2 ) 1/2 + (E\k\ 2 y/ 2 . (8.6) 

From (8.6), (7.7), (7.8) and (8.3) we obtain that 

h < C p r p + CN t ( Pt /t) 1/2 -> 0, t -> 00. (8.7) 

S'tep ii) By the triangle inequality, 

/ 2 < i|E^KI 2 -Qoo(^,^|<i|Q*(^,^)-Qoo(^,^)| 

= J 2 i + / 22 + / 23 , (8.8) 

where Q t is the quadratic form with the matrix kernel y)) . (3.5) implies that 

I21 — > . As for 7 2 2 ■ we first obtain that 

/ 22 < E \ Er l r t\- (8-9) 

\j\,\k\<N t 

The next lemma is the corollary of [13, Lemma 17.2.3]. 

Lemma 8.2 Let A,B be the subsets of Z d with the distance dist (A,B) > r > , and let 
£,?7 be random variables on the probability space (J~t a , B(H a ), A*o) ■ Let £ &e measurable with 
respect to the a -algebra a (A) , and rj with respect to the a -algebra <r(B) . Then 
1) \E£ri - E£Eri\ < Cab v? 1/2 (r) if (£|£| 2 ) 1/2 < a and (£M 2 ) 1/2 < b . 
ii) \E£r] — E£Ei]\ < Cab tp(r) if |£| < a and \r)\ < b a.s. 

We apply Lemma 8.2 to deduce that I22 — > as t — > 00 . Note that = (io(x), x^(x)$ g (-, t)) 
is measurable with respect to the a -algebra cr(it^) . The distance between the different rooms 
R\ is greater or equal to p t according to (7.1). Then (8.9) and (7.6), S3 imply by Lemma 
8.2 i) , that 

I22 < CNf^pt), (8.10) 

which tends to as t — > 00 by (8.3). Finally, it remains to check that J 23 — > , t — > 00 . 
We have 

QtVSFg, 9 g ) = E(Y , t)) 2 = E(J2M + 4) + *t) 
according to (7.5). Therefore, by the Cauchy-Schwartz inequality, 



,2. 



hs < \Eferl) -^(E^ + E^' + ^ 

< CN^EH? + Ci(^(E^') 2 ) 1/2 (^E^^| 2 + E\h\ 2 ) 1 ' 2 + ^l^l 2 . (8.11) 
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Then (7.6), (8.9) and (8.10) imply 

£(E/i) 2 < E t W| 2 + E \Efi^\ < CN t A t /t + C^y/ 2 ^) < C 2 < oo. 

j ^ J; 

Now (7.7), (7.8), (8.11) and (8.3) yield 

hs < C 1 N^p t /t + C 2 N t ( Pt /t) 1 / 2 + C 3 t- p ^0, t^oo. (8.12) 
So, all the terms I 2 i , I22 , ^23 in (8.8) tend to zero. Then (8.8) implies that 

h<l\J2E\ri\ 2 -Qoo(* g ,* 9 )\^0, t^oo. (8.13) 

Step Hi) It remains to verify that 

/3 = |i?exp^E^}-exp{-^(E^) 2 }|^0, t^oc. (8.14) 

Lemma 8.2, ii) with a = b = 1 yields: 

|Sexp{iE/t} - II Eexp{iri}\ 

-N t 

N t -1 _ N t -1 

< \E exp{ir^ Nt } exp{i E ^} — exp{?r t " Ar *}£'exp{i E r t}\ 

-N t +l -N t +1 
Nt-1 N t -1 

+ \E exp{ir^ Nt }E exp{i E r t} ~ ][ E exp{ir J t }\ 

-N t +1 -N t 
N t -1 Nt-1 

< C<p(p t ) + \Eexp{i E H}~ II Eexp{iH}\- 

-Nt+1 -Nt+1 

Then we apply Lemma 8.2, ii) recursively and get, according to Lemma 8.1, 

Nt-1 

\Eexp{i£r}} ~ II Eexp{iri}\ < CN t <p(p t ) -> 0, t -> 00. (8.15) 
It remains to check that 

JVt-l n 

I I] J Eexp^}-exp{--E^IWr}|^0, t - 00. (8.16) 

-N t 2 

According to the standard statement of the Lindeberg Central Limit Theorem (see, e.g. [16, 
Theorem 4.7]) it suffices to verify the Lindeberg condition: V5 > 



a 2 . 



Here a t = J2tE\ r t | 2 ? an d E a f := E(^X a f^j , where X a is the indicator of the event |/| > 
Note that (8.13) and (8.1) imply that a t — > Qoo(^ 9 , ^ g ) ^ 0, £ — > 00. Hence it remains to 
verify that 

~y] E a \ri\ 2 — > 0, t — > 00, for any a > 0. 

This follows from the bounds for the fourth order moments as in [6, Section 9]. This completes 
the proof of Proposition 2.12. □ 
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9 Ergodicity and mixing for the limit measures 

The limit measure fj,^ is invariant by Theorem A iv) . Let denote the integral over 

Theorem 9.1 Let all assumptions of Theorem A hold for the equation (1.1) and the initial 
measure /1q . Then U (t) is mixing with respect to the corresponding limit measure /i^ , i. e. 
V/,5 6i 2 (H a ,/ioo) 

lim E^fWYMY) = E 00 f(Y)E 00 g(Y). (9.1) 

t— >oo 

In particular, the group U (t) is ergodic with respect to the measure floo : 

1 T 

lim - / f(U(t)Y)dt = E^fiY) {mod /i^). (9.2) 



Proof Since fi^ is Gaussian, the proof of (9.1) reduces to the proof of the following con- 
vergence: V^i, ^ 2 e V 

lim E 00 (U(t)Y^ 1 )(Y^ 2 ) = 0. (9.3) 

t — >oo 

Using the Parseval identity and (4.8) we obtain similarly to (6.4) that 

£0^(^*1X^*2} = (2n)- 2d j (g(t,6) qoo (6),f (6)^,(9) ®T 2 (6))d6 

+ (2tt)- m | (g(t,6) qoo (6), 9(9)^,(9) ®T 2 (6))d6 

= If(t)+Ig(t). (9.4) 

Lemma 9.2 The uniform bound holds: \\Q(t, 0)<Zoo(0)|| < G{9) , t > 0, where G{9) E 
L\T d ) . 

Proof (4.6) implies that 

id, m,m = ( _ a SnT.n ) ( ) ■ (»■«> 

\\g{t,0)q oo (6)\\<C E (9-6) 

»J'=0,1 

It remains to prove that E L 1 (T d ) . Since ^(9) E L 1 (T d ) by Remark 2.9, it 

suffices to verify that fi~ 1 (6»)g^(6») G L 1 (T d ) , j = 0, 1 . This also follows from Remark 2.9 
if Co = . Otherwise, we will use the condition ES. Namely, owing to (2.22), we have 

^mae) = £ n^ir^M^ri^) (9.7) 

(T=l 
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since Q 1 (9) commutes with its spectral projection U a (9) . At last, (2.21) and ES imply 



fi^Mj 1 = ^(p-^ + ffiSl) G L\T d ). □ 

The Lemma 9.2 together with (3.11) and Lemma 3.3 imply that \/5 > 3e > such that 

\I f (t)\<5, t>0. (9.8) 
It remains to study the oscillatory integral I g (t) . Rewrite it using (5.1), in the form 

W = E E / gMe ±iUk{d)t a^e)(q 00 (e),if 1 (e)®¥ 2 (e)) <w. (9.9) 

m ±,fcGln Td 

Here all phase functions oj k {0) and the amplitudes aj:(0) are smooth functions in the 
suppg m . Furthermore, Vu>k(9) = only on the set of the Lebesgue measure zero. This 
follows similarly to (2.4) since VuJk{0) ^ const by the condition E4. Hence, 

I g (t) ^ as t -> oo, (9.10) 

by the Lebesgue-Riemann Theorem since G L x (T d ) . Finally, (9.4)-(9.10) imply (9.3) 
since 5 > is arbitrary. □ 
Remark Similar result for wave and Klein-Gordon equations has been proved in [4, 5]. 



10 Appendix: Crossing points 
10.1 Proof of Lemmas 2.2 and 2.3 

Step 1 By the condition El the matrix V{9) is analytic function in a connected open 
(complex) neighborhood O c (T d ) of T d in T d := T d ©iR d . Consider the analytic function 
d(9,u) := det(V(9) — uj 2 ) in O c (T d ) x C and the analytic subset defined by the equation 
d(9,u) = in O c (T d ) x C . The subset consists of the points (0,±u k (0)), k e I„ . It 
is important that d(9 ) uj) ^ for any fixed 9 G O c (T d ) , hence the function d satisfies 
the Weierstrass condition of [15, Section 2.1.1]. Therefore, by the Weierstrass Preparation 
Theorem in [15, Thm 2.1], there exists a proper analytic discriminant subset A C O c (T d ) s.t.: 
for O G O c (T d ) \ A there exists a (complex) neighborhood O c (Q) of in O c (T d ) where 
each of 0Jk{9) can be chosen as a holomorphic function. More precisely, this is established 
in the proof of [15, Proposition 2.1] which is the main step to the proof of the Weierstrass 
Theorem. We set C* := A n T d and 0(G) = C c (0) n T d for G T d \ C* . Then Lemma 
2.2 ii) follows for ujk{9) ■ 

Step 2 The identity (2.4) will follow from next general Proposition. 

Proposition 10.1 Let M. be a proper analytic subset of O c (T d ) . Then the Lebesgue mea- 
sure of the intersaction M = M. H T d is zero. 
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Proof Let us use the analytic stratification of the analytic sets which is constructed in [12, 
Thm 19 of Chapter II.E and Thm 10 of Chapter III. A]. Namely, for each G M there 
exists a complex neighborhood O c {&) s.t. M. n O c (Q) = Li <s<d~i.Ms , where each Ms is 
an analytic submanifold of the complex dimension 5 < d — 1 : here we use that M. is the 
proper analytic subset in O c (Q) . Now 

Mna(e) = u <K (l _ 1 (^nT d ). (10.1) 

Lemma 10.2 Let O G M and 5 = 0,...,d— 1. T/ien t/iere exists a (real) neighborhood 
0(Q) of 6 m T d s.t. the intersection A4$ fl 0(Q) is contained in a smooth submanifold 
of T d of the real dimension d — 1 . 

Proof We may assume that i) Ai$ is defined by the equations hj(9) — , j = 1, d — 5 , 
with the holomorphic functions foj in C (0) , and ii) V c hj(9) ^ 0, 9 G C? c (@) , where V c 
stands for the complex gradient. It is important that d — 5 > 1 so we have at least one 
function hi(9) . Then h\{9) = f\{9) + igi(9) with the real smooth functions fi,gi, and 
f\{9) = gi {9) = , 9 e M 5 f] C c (0) . However, V c h j (9) = V r + iV r gi(9) ^ , where 
V r stands for the real gradient. Therefore, either V r /i(0) ^ or V r gi(0) ^ . □ 

Now Proposition 10.1 obviously follows. □ 
This Proposition implies (2.4) since A is a proper analytic subset of O c (T d ) . Lemma 2.3 
also follows from Proposition 10.1 since E4 implies that det V{9) ^ in T d and D k {9) ^ 
in T d \ • 

Step 3 Lemma 2.2 in) follows from the construction in [15, Section 2.1]. Lemma 2.2 iv) 
follows from (2.6) since the projection U a (9) can be expressed by the Cauchy integral over 
the contour surrounding the isolated eigenvalue oj Ta (9) . 

Step 4 It remains to prove Lemma 2.2 ii) . Let O(0) denote a small real neighborhood 
of a point 0eT rf \C, and E a {9) = n CT (^)R ri . It suffices to construct an orthonormal basis 
{e k (9) : k E (r^i, 7V]} in E a (9) which depends real-analytically on 9 E 0(Q) . 

Let us choose an arbitrary basis {b k (Q) : k e (r r7 _ 1 ,r r7 ]} in E a (Q) . Then ILj-^^G) 
depend real-analytically on 9 e 0(Q) , and {n o -(6')6 / fc(6) : k e (r a -±, r a }} is a basis of 
E a {9) for 9 from a reduced neighborhood 0'{&) . Finally, construct the othonormal basis 
{e k {9) : k G (r a -±,r a ]} by the standard Hilbert-Schmidt orthogonalization process applied 
to {U a (9)b k (e) : k G (rv^rv]} for each 9 G 0'{Q) . □ 

Remark 10.3 Lemma 2.2 in) also follows from [15, Section 2.1] since the enumeration 
(2.5), (2.6) corresponds to the factorization of type [15, (2.5)] for the function d(9,u) , into 
the product of the irreducible factors, with the multiplicities r a — r„-\ , which is constructed 
in [15, Thm 2.1]. 

10.2 Proof of Lemma 2.4 

Step 1 Let us fix arbitrary k, I G I n and consider u k {9) as the functions of V G IZn and of 
9 G T d . It suffices to prove that D k (9) and V(uo k (9) ± uji(9)) are analytic and are not zero 
in an open dense subset in 1Z N x T d . 
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Let us consider Vk>i'{x) , k',l' 6 I„, < N , as the coordinates of the matrix-function 
V in the region TZn ■ Condition E2 allows us to consider Vk>i>(x) as independent real 
variables for any k',l' G l n and the points x with either x\ > , or Xi = and x<i > , 
or Xi = x 2 = and rr 3 > , etc. Let us identify TZjy with corresponding range R M of the 
independent real variables Vkn>{x) . 

Step 2 Consider u>k(9) as the functions of {V^i'^x)} and 6 in C M x T d . As above, each 
u)k{0) can be chosen as a holomorphic function outside a proper analytic discriminant subset 
A C C M x Tf . Lemma 10.2 implies that the region O : = (R M x T d ) \ A is an open dense 
subset in R M x T d . Therefore, it suffices to prove that the functions D k and V(tOk ± wi) 
are not identically zero in each connected open component of O . However, the region of 
analyticity O := (C M x Tf) \ A is connected. Hence, it remains to construct a point of 
qm x r^d g -j. functions D k and V(cjfe ± uoi) are holomorphic and non identically zero 
in a neighborhood of the point. It is easy to construct such point for any n > 1 : we can 
choose an arbitrary 9 G T d and the nearest neighbor crystal (1.12) repeated n times with 
distinct masses m k , k G I n . □ 
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